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NONLINEAR  PROBLEMS  IN  FLUID  DYNAMICS  AND 
INVERSE  SCATTERING 


ABSTRACT 

The  inverse  scattering  of  a  class  of  differential-difference  equations  and  multidimensional 
operators  has  been  constructed.  Solutions  of  nonlinear  wave  equations  arising  in  fluid  dynam¬ 
ics  and  plasma  physics  were  found  and  analyzed.  Novel  and  fast  computational  algorithms 
were  developed  and  in  relevant  nonlinear  problems  wavelet  bases  with  controlled  localization 
properties  in  the  time-frequency  domain  have  proven  to  be  an  effective  tool.  Landau  type  am¬ 
plitude  equations  governing  the  small  gap  Taylor  problem  have  been  obtained  and  analyzed. 
Detailed  project  summaries  and  research  activities  are  given  in  the  attached  document.  A 
list  of  published  papers,  preprints,  and  invited  lectures  are  included. 
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Inverse  Scattering  and  Nonlinear  Waves 
M.J.  Ablowitz 


Abstract 

Research  investigations  in  the  study  of  inverse  scattering  and  nonlinear  waves  continues 
with  a  number  of  significant  results  obtained.  Solutions  and  properties  of  physically  signifi¬ 
cant  multidimensional  nonlinear  wave  equations  and  associated  inverse  scattering  problems 
have  been  obtained  and  analyzed.  Novel  computational  methods  to  solve  nonlinear  systems 
have  been  developed  and  in  certain  cases  underlying  numerical  chaos  has  been  identified  and 
described  analytically.  An  important  new  class  of  applicable  multidimensional  nonlinear 
equations  whose  solutions  become  unbounded  in  finite  time  has  been  analyzed. 

Overview 

The  past  year  has  been  a  very  active  period  for  our  research  program.  A  total  of  eight 
papers  have  been  published,  one  book  which  summarizes  the  recent  developments  in  the 
field  has  been  printed  by  Cambridge  University  Press  and  eight  preprints  have  either  been 
accepted  or  have  been  submitted  for  publication.  Our  preprints  are  also  issued  as  reports  of 
the  Program  in  Applied  Mathematics  at  the  University  of  Colorado  before  they  are  published. 
An  overview  and  brief  description  of  our  research  investigations  is  given  below.  Complete 
details  may  be  found  in  our  research  papers.  A  list  of  recent  reprints  and  preprints  is 
compiled. 

Research  Activities 

Multidimensional  Nonlinear  Wave  Equaitons  and  Inverse  Scattering 

Our  investigations  of  physically  interesting  nonlinear  wave  equations  continues.  The  deep 
connection  between  integrable  nonlinear  wave  equations  and  inverse  scattering  allows  us  to 
analyze  and  solve  nonlinear  systems  which  otherwise  would  be  intractable.  One  example 
of  a  particularly  natural  and  significant  system  is  the  two  space-one  time  extension  of  the 
well-known  Toda  lattice.  Using  inverse  scatteimg  methods  we  have  been  able  to  solve  this 
system  in  a  variety  of  cases.  In  the  case  where  the  linearized  system  is  analogous  to  an 
elliptic  ODE  we  can  formulate  and  solve  the  relevant  boundary  value  problem.  In  one  case 
we  have  discovered  a  nonlinear  analgoue  of  the  Sommerfeld  radiation  condition  (frequently 
required  in  linear  problems)  which  must  be  imposed  in  order  to  obtain  a  unique  solution.  It 
is  necesary  to  impose  such  a  radiation  condition  in  a  much  wider  class  of  problems. 

Reductions  of  certain  physically  interesting  multidimensional  systems  has  led  us  to  study 
a  class  of  nonlinear  ODE’s  which  have  novel  features.  The  prototypical  situation  is  character¬ 
ized  by  a  system  of  three  coupled  ODE’s  first  discovered  by  Halphen  in  1881.  We  transform 
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this  system  to  a  third  order  sealar  equation  studied  by  Chazy  in  1910.  We  have  shown  that 
Chazy’s  classical  solution  can  in  fact  be  expressed  in  terms  of  autornorphic  functions,  and  in 
particular  well-known  modular  forms.  We  have  also  found  generalizations  of  this  system  for 
which  the  solution  can  be  expressed  in  terms  of  more  general  autommorphic  functions.  As 
such  these  ODE’s  provide  an  important  new  class  of  “  special  functions”  arising  in  integrable 
systems  which  are  sure  to  be  heavily  studied  by  researchers  in  this  field. 

Computational  Algorithms  Instabilities  and  Numerical  Chaos 

We  continue  our  studies  which  have  already  led  to  a  number  of  valuable  discoveries.  We 
have  shown  that  different  numerical  discretizations  of  integrable  nonlinear  wave  equations 
can  lead  to  dramatically  different  results,  even  though  the  discretizations  differ  slightly.  For 
example  modifications  of  how  one  discretizes  the  nonlinear  terms  in  the  system  can  lead  to 
joint  spatial  and  temporal  numerical  chaos.  Based  upon  extensive  numerical  experiments 
we  have  identified  two  excellent  numerical  schemes.  One  of  them  is  a  differential-difference 
equation  suggested  by  inverse  scattering  analysis.  The  other  is  the  well-known  Fourier  split 
step  algorithm.  The  differential-difference  algorithm  is  easy  to  implement,  is  effective  and 
is  considerably  faster  than  Fourier  algorithms.  Nevertheless,  even  with  the  “best”  schemes 
available  we  find  that  associated  with  a  typical  class  of  initial  conditions,  temporal  chaos 
can  be  excited  by  roundoff  errors.  In  fact  the  underlying  solutions  are  linearly  unstable  in 
this  regime  of  parameter  space  and  this  causes  miniscule  errors  to  grow  rapidly.  The  errors 
eventually  become  as  large  as  the  main  wave  even  though  the  schemes  preserve  the  conserved 
quantities  quite  well.  These  ideas  are  being  examined  in  the  context  of  many  other  physically 
interesting  nonlinear  equations. 

Focusing  Singularities  in  Multidimensional  Nonlinear  Wave  Equations 

We  have  demonstrated  that  a  certain  class  of  multidimensional  extensions  of  the  well- 
known  Korteweg-de Vries  equations,  often  referred  to  as  higher  nonlinear  Kadomtsev- 
Petviashvili  (KP)  equations,  have  solutions  which  will  become  infinite  in  finite  time.  This 
phenomenon  is  sometimes  referred  to  as  “wave  collapse.”  We  have  examined  the  behavior  of 
the  generalized  KP  system  both  analytically  and  numerically,  and  have  shown  that  above  a 
certain  threshold  for  the  power  of  the  relevant  nonlinear  term  there  can  be  a  focusing  singu¬ 
larity;  e.g.  this  occurs  for  the  so  called  modified  KP  equation  which  arises  in  plasma  physics. 
This  system  is  yet  another  well-known  nonlinear  wave  equation  which  possesses  blow  up 
singularities.  Other  multidimensional  nonlinear  wave  equations  which  have  this  property 
include  the  two  dimensional  nonlinear  Schrodinger  equation  and  the  Davey-Stewartson  sys¬ 
tem.  Each  of  these  equations  have  been  extensively  studied  by  many  researchers  in  a  variety 
of  fields  of  application  including  nonlinear  optics,  fluid  dynamics,  and  plasma  physics. 
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Propagation  and  capturing  of  singularities  in  problems  of  fluid 

dynamics  and  inverse  scattering 
Gregory  Beylkin 


Abstract 

The  study  and  development  of  the  fast  adaptive  algorithms  for  solving  nonlinear  equa¬ 
tions  using  orthonormal  bases  with  controled  localization  in  the  time-frequency  domain  (e.g. 
wavelets)  has  been  succesftil  in  establishing  an  approach  for  the  use  of  the  wavelet  bases 
in  numerical  computations.  This  approach  is  novel  and  involves  the  use  of  sparse  exponen¬ 
tials  for  the  evolution  in  time  as  well  as  two  key  algorithms  for  adaptive  differentiation  and 
multiplication  in  the  wavelets  bases.  Several  additional  algorithms  relevant  to  solving  the 
nonlinear  equations  and  capturing  the  singularilies  (shocks)  are  at  various  stages  of  develop¬ 
ment.  An  overview  and  a  brief  description  of  the  research  effort  is  given  below.  The  details 
may  be  found  in  the  research  reprints  and  preprints  compiled  for  this  report. 

Research  Activities 

During  this  year  the  research  and  development  of  numerical  algorithms  for  solving  non¬ 
linear  equations  grew  into  a  sizable  effort.  I  work  now  with  two  graduate  students,  .lames 
Reiser  and  Robert  Cramer.  With  Robert  Cramer  we  are  working  on  the  multidimensional 
algorithms  for  the  fast  evaluation  of  a  class  of  integral  operators  on  functions.  Such  algo¬ 
rithm  will  be  useful  in  evaluating  integral  operators  which  appear  in  fluid  dynamics  (e.g. 
Biot-Savart  law  relating  vorticity  and  velocity  of  the  fluid).  With  James  Reiser  we  are 
developing  adaptive  algorithms  for  the  evolution  of  nonlinear  equations.  So  far  we  have 
worked  on  equations  in  one  spatial  dimension  in  order  to  test  our  approach  and  evaluate  its 
advantages/disadvantages  in  comparison  with  other  methods.  We  pay  special  attention  to 
ensure  that  the  algorithms  we  are  developing  generalize  to  multiple  dimensions  in  an  effec¬ 
tive  manner.  In  addition  to  this  effort,  several  algorithms  were  developed  (jointly  with  M.E. 
Brewster)  which  use  wavelet  bases  on  the  interval,  thus  allowing  us  to  incorporate  boundary 
conditions  properly. 

Numerical  Solution  of  Nonlinear  Equations  in  Wavelet  Bases 

On  one  hand,  solving  partial  differential  equations  in  wavelet  bases  has  some  features 
similar  to  those  of  pseudo-spectral  methods,  multigrid  methods  and  adaptive  grid  methods. 
On  the  other  hand,  solving  in  wavelet  bases  has  many  new  features  that  permit  fast  and 
accurate  numerical  algorithms.  We  have  developed  several  algorithms  which  allow  us  to 
solve  a  class  of  nonlinear  partial  differential  equations. 

A  fast  adaptive  algorithm  for  computing  the  pointwise  product  of  functions  represented 
in  wavelet  bases  has  been  developed.  The  number  of  operations  is  proportional  to  the 
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complexity  of  the  representation  of  the  function  in  a  wavelet  basis.  e.j>.  for  smooth  functions 
with  a  finite  number  of  singularities  the  algorithm  requires  (V(log.Y)  operations,  where  log  .Y 
is  the  number  of  dyadic  scales.  The  constants  in  tin*  complexity  estimate  are  reasonable  and 
the  algorithm  is  easily  extendable  to  the  multidimensional  ease. 

A  fast  adaptive  algorithm  for  computing  the  derivatives  in  wavelet  bases  using  the  non¬ 
standard  form  has  been  developed.  This  algorithm  is  applicable  to  a  class  of  operators  whi  -h 
includes  some  singular  integral  operators  such  as  the  Hilbert  Transform. 

For  evolution  in  time,  we  compute  the  exponential  of  an  operator  (e.g..  of  the  second 
derivative  operator)  in  wavelet  bases.  Such  exponentials  are  sparse  in  wavelet  bases  and 
result  in  an  unconditionally  stable  schemes.  The  accuracy  of  this  approach  is  superior  to 
other  time  stepping  schemes  since  those  attempt  to  approximate  the  exponential. 

So  far  we  used  these  algorithms  to  solve  Burger’s  equation.  We  plan  to  apply  this 
approach  to  several  nonlinear  equations  in  order  to  evaluate  the  performance  of  our  method. 

Fast  Algorithms  for  Evaluation  of  Multidimensional  Integral  Operators 

We  are  implementing  a  fast  algorithm  for  evaluation  of  a  class  of  multidimensional  integral 
operators  on  functions,  where  special  care  has  been  taken  to  assure  small  constants  in  the 
complexity  estimates.  We  expect  to  compare  the  speed  of  the  new  algorithm  with  that  of 
the  Fast  Multipole  Method. 

Two-point  Boundary  Value  Problems 

A  hybrid  method  for  solving  a  two-point  boundary  value  problem  has  been  developed. 
The  problem  is  formulated  using  ordinary  finite  difference  approximation  and  then  is  solved 
in  a  wavelet  basis.  The  main  tool  in  this  approach  is  the  diagonal  preconditioning  available 
for  the  periodized  differential  operator  in  the  wavelet  bases.  As  a  result,  the  fundamental 
solution  of  the  two-point  boundary  value  problem  for  the  elliptic  differential  operators  has 
been  computed  in  O(N)  operations. 

Availability  of  the  fundamental  solution  as  a  sparse  O(N)  matrix  allows  us  to  take  a  dif¬ 
ferent  approach  to  the  implicit  schemes  for  solving  evolution  equations  (e.g.  Crank-Nicolson) 
by  inverting  the  operator  (matrix)  and,  thus,  making  them  explicit.  We  are  looking  at  several 
applications  to  take  advantage  of  such  a  procedure. 

Wavelet  Bases  on  the  Interval,  Multiresolution  Numerical  Algorithms  and  Ho¬ 
mogenization 

We.  are  starting  to  apply  the  wavelet  bases  on  the  interval  as  a  means  t,o  solve  the 
boundary  value  problems.  We  expect  to  construct  high-order  multiresolution  methods  for 
evolution  in  time.  We  are  also  looking  into  the  extension  of  functions  from  the  interval 
and  splitting/patching  of  functions  defined  on  an  interval  with  possible  signal  processing 
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applications. 

Joinrlv  with  Mary  E.  Brewster,  we  have  constructed  a  numerical  scheme  foi  homoge¬ 
nization  of  linear  differential  and  difference  equations  via  mult  i  resol  tit  ion  analysis.  Though 
this  work  is  in  its  beginning,  it  has  a  potential  for  a  significant  impact  by  permitting  us 
to  construct  the  equations  for  computing  the  projection  on  a  sparse  scale  of  the  solution  of 
equations  given  on  the  very  fine  scales.  Om  homogenization  method  may  be  used  as  part 
of  a  multiresolution  strategy  to  obtain  the  Green's  function.  In  this  approach  the  Green's 
function  is  constructed  on  successively  finer  scales.  Combining  the  multiresolutiou  strategy 
with  wavelet  bases  on  the  interval  results  in  a  new  class  of  high-order  methods  for  evolution 
equations. 
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Langmuir  Circulation  and  Spiral  Flows 
Duane  P.  Sather 


I.  Summary  of  Research  Activities. 

(1)  Landau  type  amplitude  equations  for  the  small-gap  Taylor  problem  were  derived 
and  analyzed  in  [2]  (see  the  list  of  publications  in  II)  The  equation  obtained  are  global 
and  are  more  complete  than  those  determined  by  a  formal  two-timing  analysis.  Transition 
solutions  (i.e.,  heteroclinic  orbits)  connecting  the  trivial  Couette  flow  with  bifurcating 
steady  flows  were  obtained  by  solving  singular  evolution  equations  in  minute- dimensional 
spaces.  The  results  obtained  lead  to  the  first  analytic  description  of  the  difference  between 
primary  and  secondary  flows  of  a  viscous  fluid  and  provide  a  more  detailed  description  of 
the  transition  solutions  than  that  obtained  by  the  use  of  formal  two-timing  methods  or 
center  manifold  theory. 

The  work  in  (2)-(4)  is  joint  with  George  H.  Knightly  of  the  University  of  Mas¬ 
sachusetts. 

(2)  The  bifurcation  and  stability  properties  of  spiral  flows  for  rotating  plane  Couette 
flow  problems  were  obtained.  These  are  the  first  complete  results  on  the  bifurcation  and 
stability  of  viscous  spiral  flows.  The  final  version  of  these  results  is  published  in  [3] . 

(3)  Periodic  waves  were  shown  to  exist  and  their  expansions  obtained  in  [1;4]  for 
various  Couefcte-Poiseuille  problems  for  viscous  spiral  flows.  The  existence  of  periodic 
waves  in  such  problems  has  been  conjectured  but  never  proved. 

(4)  A  continuum  of  periodic  waves  bifurcating  supercrifcically  from  the  basic  spiral  flow 
was  obtained  in  [5]  for  rotating  plane  Couette  flow.  Results  of  this  type  provide  a  simple 
explanation  for  the  occurrence  of  turbulent-like  flows  in  spiral  flow  problems  and  also  the 
first  analytic  description  of  such  turbulent-like  flows.  The  usual  Hopf  bifurcation  theory 
does  not  apply  to  such  problems  and  new  perturbation  methods  using  singular  operators 
were  developed. 

(5)  The  basic  equations  for  Langmuir  circulations  in  upper-ocean  mixing  problems 
when  the  Stokes  drift  has  a  cross-wind  component  were  developed  in  [6];  the  set  of  partial 
differential  equations  are  derived  from  a  rational  model  that  describes  the  development  of 
mean  currents  and  Langmuir  circulations  as  a  single  system  driven  by  a  prescribed  wind 
stress  and  a  prescribed  surface  wave  field.  In  physical  situations  where  there  is  no  cross- 
wind  component  of  the  Stokes  drift  and  one  seeks  only  roll- like  solutions  independent 
of  the  wind  direction,  the  equations  in  [6]  reduce  to  those  derived  previously  by  other 
investigators.  The  equations  in  [6]  are  derived  by  the  use  of  multiple  time  scales  related 
to  the  interaction  of  surface  waves,  wind  and  currents.  The  results  in  [6]  provide  the  only 
rational  theoretical  explanation  to  date  of  observed  phenomena  in  Langmuir  circulations 
such  as  cross-drift  of  Langmuir  cells  and  the  rapid  breakdown  of  Langmuir  cells  under  a 
change  in  wind  direction. 
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(6)  A  mathematical  analysis  of  the  equations  derived  in  [6]  is  developed  in  [7].  The 
usual  methods  of  Hopf  bifurcation  provide  only  partial  results  here  and  new  perturbation 
methods  related  to  those  in  (4)  using  singular  operators  were  developed  to  treat  Langmuir 
circulation  problems.  Such  methods  provide  a  more  detailed  description  of  the  resultant 
periodic  waves  than  those  obtained  by  the  use  of  formal  two-timing  methods  or  center 
manifold  theory.  Other  problems  for  Langmuir  circulations  involving  phenomena  such  as 
changes  in  the  number  of  cells  and  patterns  other  than  rolls  also  can  be  treated  using  tie- 
basic  model  equations  derived  in  [6].  Such  problems  will  be  considered  in  a  series  of  papers 
in  preparation. 

II.  Publications  and  Papers  in  Preparation. 

The  following  papers  were  completed  and  accepted  for  publication  during  the  period 
of  the  Progress  Report.  Preprints  and  reprints  have  been  forwarded  to  ONR. 

[1]  “Time-periodic  states  in  problems  containing  a  structure  parameter,”  (with  George  H. 
Knightly),  Progress  in  Partial  Differential  Equations  (J.  Bandle,  et  al.,  eds.),  Pitman 
Research  Nvtes  in  Mathematics  Series,  John  Wiley  &  Sons,  Inc.,  New  York,  1992,  pp. 
34-92. 

[2]  “Transition  solutions  in  the  Taylor  problem,”  Arch.  Rational  Mech.  Anal.  121  (1992), 
267-301. 

[3]  “Periodic  waves  in  rotating  plane  Couette  flow,”  (with  George  H.  Knightly),  Z.  Angew. 
Math.  Phys.  44,  1-16. 

[4]  “Symmetry  in  rotating  plane  Couette-Poiseuille  flow,”  (with  George  H.  Knightly), 
Lectures  in  Applied  Mathematics,  American  Mathematical  Society,  accepted  for  pub¬ 
lication,  refereed  article,  11  pages. 

The  following  papers  were  completed  and  submitted  for  publication  during  the  period 
of  the  Progress  Report.  Preprints  have  been  forwarded  to  ONR. 

[5]  “Continuaof  periodic  waves  in  rotating  plane  Couette  flow,”  (with  George  H.  Knightly), 
J.  Diff.  Eqns.,  submitted  for  publication,  34  page  manuscript. 

[6]  “Langmuir  circulations  when  the  Stokes  drift  has  a  cross-wind  component,”  (with 
George  H.  Knightly),  J.  Fluid  Mech.,  submitted  for  publication,  26  page  manuscript. 

The  following  paper  is  in  the  final  stages  of  preparation.  A  preprint  will  be  forwarded 
to  ONR  when  it  is  available. 

[7]  “Periodic  traveling  waves  in  Langmuir  circulations,”  (with  George  H.  Knightly),  to 
be  submitted  to  Arch.  Rational.  Mech.  Anal. 
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A  *  denotes  invitations  accepted 

(1)  Applied  Dynamics  and  Bifurcation  Theory,  Oberwolfach,  Germany,  January  12  IS. 

1992. 

(2)  First  World  Congress  of  Nonlinear  Analysis,  Florida  Inst,  of  Technology,  Melbourne, 
August  19-26,  1992. 

(3) *  Summer  Seminar  in  Applied  Mathematics,  “Exploiting  Symmetry  in  Applied  and 

Numerical  Analysis,”  sponsored  by  AMS-SIAM,  Fort  Collins,  July  26- August  1,  1992. 

(4) *  Chapman  Conference  on  “Double-Diffusive  Convection,”  sponsored  by  the  American 

Geophysical  Union,  Phoenix,  November  3-6,  1993. 


13 


